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^ ^ Abstract. We study the limiting behavior of viscous incompressible 

flows when the fluid domain is allowed to expand as the viscosity van- 
ishes. We describe precise conditions under which the limiting flow 
, satisfles the full space Euler equations. The argument is based on trun- 

cation and on energy estimates, following the structure of the proof of 
Kato's criterion for the vanishing viscosity limit. This work comple- 
, ments previous work by the authors, see |21 |8j . 
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1. Introduction 

In [5], the second and third authors, in collaboration with Drago§ Iftimie, 
showed that, if an obstacle is scaled by a factor e, then in the limit as 
viscosity vanishes the solutions to the Navier-Stokes equations external to 
the obstacle converge strongly in L°°([0, T]; L^) to a solution to the Euler 
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equations in the whole space, as long as e < az/ for a specific constant a. 
They also give the rate of convergence in terms of v and e. 

In [8], the first author considered the complementary problem of large 
domain asymptotics, studying convergence to full plane flow of solutions of 
Euler or Navier-Stokes in a large domain. The present article is a natural 
continuation of both [5] and [8] . 

For a domain with boundary, it a classical open problem whether solutions 
of the Navier-Stokes equations converge to solutions of the Euler equations 
when viscosity vanishes. In [5j the authors are considering two limits si- 
multaneously: the vanishing viscosity limit and the limit as the obstacle 
shrinks to a point, solving the external problem for the Navier-Stokes equa- 
tions. This means studying the way in which a small boundary obstructs 
the vanishing viscosity convergence. Here, we consider what happens as a 
bounded domain expands by a factor R to fill the whole space, giving the 
convergence rate in the vanishing viscosity limit for the internal problem in 
terms of v and R. In the same spirit as [5], the present work regards the 
effect of distant boundaries in the vanishing viscosity limit. 

More precisely, let J7 be a simply connected bounded domain in M^, d = 2 
or 3, with C^-boundary T and let il/j = RQ. and = RT = OUr, where 
we assume that the origin lies inside O. 

A classical solution (u,p) to the Euler equations without forcing in all of 
satisfies 

{dtu + u-Vu + \/p = in (0, T) x R*^, 
divn = in [0,r] x R*^, 

u = uo on {0} X R*^, 

where divn'' = 0. A classical solution {u'^'^,p'^'^) to the Navier-Stokes 
equations without forcing on Qji satisfies 

dtu"^^ + n*^'^ • Vn^'^ + Vp"^^ = uAu"^^ in (0, T) x Or, 

divn'"'^ = m[0,T]xnR, 

u"^^ = on [0, r] X Fr, 

u^^R = on {0} X VLr, 

where u^^ = on Fr. 

We will work, however, with weak solutions to the Navier-Stokes equations 
(to avoid having to deal with the dependence of the time of existence of 
solutions on viscosity). 

We consider the classical functions spaces, 

V{^pi) = {n G H^{^ii): divn = in ^Ir and n = on Fr| , 

H^flji) = |n G L?{^ii) : divn = in Qr and n • n = on Fr| , 

where n is the outward directed unit normal vector field to Fr. 
The spaces V(U.'^) and H{W^) are analogously defined. 
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We define the space 

VciM."^) = {n G y(R'^) : supp( curl u) is compact} . 

We will use the notation uj{u) = curl u for the vorticity associated to a given 
velocity u. 

In dimension two the condition u € Vc'(R^) requires that the total mass 
of the vorticity be zero, see Section 3.1.3 of [13] for a discussion. Hence, 
if we want to allow vorticities with distinguished sign we must allow for 
infinite energy. To this end we recall the affine spaces Em, introduced by 
J.-Y. Chemin in [1], following a construction by R. DiPerna and A. Majda, 
see [2]. We say that u £ Em li u = v + a for some v G H{^) and for 
some stationary solution a of the Euler equations whose vorticity is smooth, 
compactly supported and has integral m. More precisely, for a given by 

a = a{x) = -. — 12 / sip{s)ds, (1-1) 
fI Jo 

for some (p G C^(M+) and 27r J ip{s)sds = m. Given the arbitrariness 
in the choice of (p we will assume, without loss of generality, that (p is of 
distinguished sign. Above we used the notation x"*" = {—X2,xi) if a; = 
{xi,X2)- Notice that a;((T)(x) = v^d^^l)- 

The classical well-posedness results for weak solutions to (E) for uq in 
subspaces of £"0 = -ff(M^) remain true when Eq is replaced by Em', see, 
for instance. Theorem 5.1.1 p. 85 of [Ij (Yudovich's theorem). In particu- 
lar, fixing a value of T > 0, if uq is in Em H H^{M.'^) with compactly sup- 
ported initial vorticity then the solution u to (E) will lie in C{[0,T]; Em) H 
L°°([0,r]; 7^1(1^2)). 

Throughout this paper we will assume that the initial velocity uq for 
solutions to (E) lies in C'^(M'^) for s > 1 so that a unique solution u to the 
Euler equations {E) with initial velocity uq exists in the space C*([0, T] x 
M.'^) for all T < T*; see, for instance, Theorem 4.2.1 p. 77 of ^ (or see 
Theorem 17. H below). The time T* can be assumed to be arbitrary in two 
dimensions, where we also assume that uq lies in Em (see Theorem 4.2.4 p. 
82 of [Ij), but only finite time existence is known in three dimensions. We 
assume that the initial vorticity is compactly supported with its support 
contained in a ball of radius Rq and define 

R{T) = inf {r: suppu;(n) C [0,r] x 5,.(0)} . 

r>0 

That R{T) is finite follows in two dimensions from the transport of vorticity 
by the flow associated to u, u being bounded uniformly over finite time. But 
R(T) is also finite in three dimension, as we show in Theorem 17.11 

Definition 1.1 (Classes of initial velocities). Let s > 1. We treat the 
following three classes of initial velocities: 

I. uo is in C'{R'^)nVc{R'^), 
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II. uq is in C (M?) n EmH {M.'^) , the support of ti;(no) is compact, and 

Qi is a disk, 
III. Uq is in C«(M3) n VciM.^). 

We assume that the initial velocity Uq^ is in For such initial 

velocities it is a classical result of Leray that there exists a weak solution 
u'^'^ to the Navier-Stokes equations (NS); in two dimensions this solution is 
unique, a result due to Ladyzhenskaya. In three dimensions, global-in-time 
existence is known, but not uniqueness, so we arbitrarily choose one such 
solution for each value of v. 

Our main result is the following: 

Theorem 1.2. Let uq be in one of the three classes of initial velocities in 
Definition \l.l\ and set F{y,R) = \\u^^ — uoWL'^inn)- -^^r all T < T* there 
exists a constant C = C{s,T,0,,uq) > such that 

(1) ifs>l, 

h"^'' - ^IIlo«([o,t];L2(c«)) < (Ciu^/'' + R-^) + F{u, R)) e^^; 

(2) ifs>2, 

- n||ioo([o,T];L2(n«)) < {C{u + + F{u, R)) e^^, 

for all sufficiently large R. 

The exponent a is defined for each of the three cases as follows: 
I. a = 1, 
II. a = 1/3, 
III. a = 1/2. 

v Ft 

Of particular interest is when we define Mq' , independently of v, to be 
that unique divergence-free vector field tangent to the boundary of ^Ir whose 
vorticity on Qji is the same as that of uq. In Section [10] such a vector field 
is denoted = WrUq. We will see in Corollarv 110.21 that. 

H'"" - uoh^inn) = R) ^ HR) < (1-2) 
for all R > 2Rq, with a defined as in Theorem 11.21 In this case, the term 
F(R) in the bounds in Theorem [L2] is dominated by the other term and so, 
in effect, it disappears. 

It follows immediately from Theorem 1 1 . 2 1 that . as long as R = R{v) oo 
as 1/ ^ and F{u, R) ^ as R oo, - u\\loo(_[q^t].l2(^q^^^^^) ^ as 

1/ ^ 0. 

It was shown in [8] for Case I that if is the solution to the Euler 
equations on ^r with initial velocity Truq then Ijii'^ — m|| j^oo([o,T];L2(Cfl(^))) ~^ 
as R ^ oo. Here, Tr is a truncation operator, which will be defined 
precisely in Section [3l see ()3.2p . This result is extended in [9] to cover Case 
II and to use the projector Pv{nii) — restriction to Qr followed by projection 
into V{Qr) — in place of Tr. This gives the following corollary: 
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Corollary 1.3. Let T < T* and set Uq = Py(^q^-^uq. Then, for uq as in 
Cases I or II, 

as long as R = R{y) ^ oo as i' ^ 0. 

The energy argument in our proof of Theorem 11.21 fohows fairly closely 
the argument in [S], which itself is closely connected to Kato's argument in 
[^. We can describe in a unified way the approach of all three papers — [7|, 
[5], and this one — as follows. Let u^^ be the solution to (NS) in a domain 

and let be the solution to {E) either in the whole space or, as in [7], in 
Q itself. In [7j, O is a fixed bounded domain; in is an external domain 

which is scaled to a point by a parameter e; for us, is a bounded domain 
which is scaled by a parameter R to fill the whole space. 

Define a correction velocity u*^ to such that u'~' = on and is 
equal to zero outside a boundary layer Fs of width 6. In [7|, (5 = Cu; in 
p], (5 = e; in this paper, 6 = CR°^. Let = — vP be an "approximate 
solution" to (-E), and observe that ti"^ = on 

The goal is to bound the norm of u^^ — = u^^ — w^ — u'~" in the space 
X = L°°([0,r];L2(Q)). To do so, one first shows that \\vP\\x ^ as (5 ^ 
or oo as the case may be. Then one bounds W = u^^ —u^ in X by making an 
energy argument, the nature of the argument differing in each case. Because 
ti"^ = on no troublesome boundary terms appear, though certain other 
terms appear because is only an approximate solution to {E). 

Kato's energy argument in |7j is designed to estimate all of the uncon- 
trollable terms by the quantity 

l|Vn-1lL(re.)' (1-3) 

which, by the most basic energy argument for solutions to (NS), is bounded 
uniformly for all T and must vanish if the vanishing viscosity limit is to hold. 
Kato's innovation is to show that the vanishing of this term is sufficient for 
the vanishing viscosity limit to hold. 

The results achieved in the three papers differ most fundamentally be- 
cause for Kato Vu'-" scales like 1/z^, which is detrimental (but unavoidable), 
introducing terms into the energy argument that cannot quite be controlled. 
For us, Vn*-^ scales like which allows us to control all of these terms. 
In [5], Vn*-" scales like 1/e, but the domain shrinks in area like e^, which 
largely counteracts the detrimental effects of Vti*^. 

The research presented here is part of a series of papers aimed at studying 
asymptotic behavior of incompressible flows under singular domain pertur- 
bations. The first result in this line of research concerned ideal 2D flow in 
the exterior of a small obstacle, see [3], followed by a study of viscous 2D 
flow in the same limit, see [4]. Beyond these, this research has included ideal 
2D flows in bounded domains with multiple holes one of which vanishes, see 
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|12j . ideal or viscous 2D flow in a large domain, see [8j, 3D viscous flow in 
the exterior of a small obstacle, [6] and, most recently, 2D flow exterior to a 
smooth obstacle approaching a segment of a curve, see |10j for the ideal flow 
case and pJJ for the viscous case. The classical open problem of vanishing 
viscosity in the presence of boundaries motivated the coupling of singularly 
perturbed domain problems with vanishing viscosity, speciflcally when the 
boundary disappears as viscosity vanishes. The flrst result in this direction 
was obtained in [5j for the small obstacle limit and the current work can be 
regarded as a natural continuation of [8] in the same spirit. 

This paper is organized as follows: Section 2 contains certain notation 
we use and conventions we follow. In Section [3] we describe an approximate 
solution to the Euler equations on il/j which we use in Section [5] to 
prove Theorem 11.21 The proof of Theorem 11.21 relies, however, on a long 
series of estimates involving u^, which require us to understand how to 
take a divergence-free vector fleld deflned in the whole plane or space and 
"truncate" it in such a way that it is unchanged in the central part of the 
domain fi^j, vanishes on the boundary of Or, and yet differs in the pertinent 
norms on as little as possible from the original vector fleld. We describe 
the two dimensional version of such a truncation operator in Section [5] and 
use it in Section [6] to deflne and obtain the necessary estimates on u^. 

The deflnition and analysis of the truncation operator in three dimensions 
is markedly different from that in two dimensions. In Section [7] we derive 
uniform-in-time bounds on the decay of the velocity and its gradient for a 
solution to {E). We then deflne the truncation operator in three dimensions 
in Section [8] and obtain the estimates on in three dimensions in Section [H 
In Section [10] we prove (jl.2p . In Section [TT] we make some comments and 
state a couple of open problems. 



The symbol C stands for a positive constant that can hold different values 
on either side of an inequality, though always has the same value on each 
side of an equality. 

For a scalar function / in two dimensions we write V^/ := (— ^2/, 
In two dimensions we deflne the vorticity of a vector fleld u to be the scalar 
curl, Lo = Lo{u) := div? — d2U^ = V"*- • u. In three dimensions, we deflne the 
vorticity to be a; = lo{u) := curln; that is, lo is the three- vector. 



It is sometimes convenient in three dimensions to view the vorticity as 
the anti-symmetric 3x3 matrix A = A(n) whose entry in the i-th row, k-th 
column is ujI = ujI{u) := {dkU^ — diU^)/2. Thus, 
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Observe that the L^-norms of A and lv are equivalent, differing only by a 
multiplicative constant. 

Given a divergence- free vector field u on let u) = u}{u) be its vor- 
ticity, which we assume to have compact support. We define the associated 
two-dimensional stream function -0 as 

tP = tl;(x) = ^l \og\x -y\uj{y)dy, (2.1) 

so that Aijj = Lo and u = SJ-^il). 

Given a divergence- free vector field u on M'^ with compactly supported 
vorticity uj = u){u) we define the associated three-dimensional (vector- valued) 
stream function ^ as 

m = m{x) = / ^^c^(y) dy. (2.2) 
4vr JiRS \x - y\ 

Hence, — A^* = uj and u = curl ^ , the latter statement following since 
div ^ = 0, which in turn can be seen from the equation Adiv ^' = 0, 
div ^ ^ at oo. 

We note in passing that an alternative to this vector-valued stream func- 
tion is to define the matrix-valued stream function 

2vr 7k3 \x - y\ 

which has the property that = '}2k (^kipik- The advantage of defining the 
stream function in this way is that it can be generalized to higher dimensions. 

3. Approximate solution to the Euler equations 

Define a cutoff function ip^ in two dimensions as follows. Fix 9 in [0,1]. 

(We will ultimately choose a value of 6 that optimizes the convergence rate 

in Theorem 11.21 ) Let 5i = 1/2k, where k is the maximum curvature of 

r = di}. Let S/j be a tubular neighborhood of L/j in fi/j of uniform width 

SiR^ for all R in [I, 00). (Decrease the value of Si if necessary to insure that 

the origin is not contained in S/j.) Put coordinates (s,r) on S/j, where s is 

arc length along F, which locates a point on F, and r is the distance along 

the inward normal at that point. 

Let g in C°°{[0, Si]) taking values in [0, 1] be defined so that g{0) = g'{0) = 
"A 1 

^ ' . Then define ip^ in C°°(Or) by ip^{s, r) = g{R-'^r) 



and g = I on 



for points (s, r) in S/j, and cp^ = 1 on $7/? \ S/j. Observe that 

IML^in,) < CR^', II VVv.^Loo(f.,) < CR-''^ (3.1) 

and similarly for higher derivatives of ip^, where C is independent of R in 
[1, 00), and tp^ = ^ and V^p^ = on F. 

We define tp^ in three dimensions more simply. Let 

S = {x G O : dist(x,F) < 1/2k} , 
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where k is the maximum of all sectional curvatures over all points of T. Let 
ip in C°°(r2) taking values in [0, 1] be defined so that = 1 on 17 \ S and 
= 0, V(/9 = on r, and let ip^{-) = if{-/R) and T,r = RT.. Then (lOl 
holds with 6 = 1. 

Let ip be the two-dimensional stream function associated to the full-plane 
Euler velocity u, as in ()2.ip . We define the vector field on O/j by 

= Tru := V^iif^ip). (3.2) 

Notice that this defines an operator Tr whose properties we will explore 
later. 

If ^ is the three-dimensional stream function associated to the full-space 
Euler velocity n, as in (|2.2p . then we define the approximation on Qr by 

= Tru := V X {ip^'^). (3.3) 

The operator Tr in both cases has the property that = Tru lies not 
just in H{Qr) but in V{Qr), and so vanishes on the boundary. It also 
satisfies (E) in 0,r \ T,r. In this sense, it is an approximate solution to {E). 

Clearly, satisfies the identity 

dtu^ = -Lp^u ■ Vu - <p^Vp + dtipW-^Lp^ (3.4) 

in two dimensions and 

dtu'^ = -(p^u-Vu- (p^Vp + Wif^ X dt^, (3.5) 

in three dimensions. 

Next we state a proposition which contains the key estimates on that 
we will use in Section[4]to prove Theorem II. 21 We prove the two-dimensional 
case of this proposition in Section [6] and the three-dimensional case in Sec- 
tion [3 

Proposition 3.1. For all T < T* , for all sufficiently large R, we have 

(1) ||Vn^||^oo([o,T];L2{Qii)) < C, 

(2) \\u^\\l°°([q,t]xQ.r) < C, 

(3) ||V?x^||ic»([o,T]xf7B) ^ ^' 

(4a) ||pV(/?-^||ioo([o,T];L2(nfl)) + l!'9tV'V99^||ioo([o,T];L2(nfl)) 

< CR-^ in 2D, 

(4b) \W<P^L'^{[Q^T];L^{nR)) + l|Vl/3-^ X ^ H^oo ([o,T];L2(Qfl)) 

< CR-^ in 3D, 

(5) ||An^||icx=([o,r];L2(nfl)) < C when s>2, 

(6) ||n^ - n||ioc([o,T];L2(Q^)) + \\u^ - y^^u\\Loo(^[o^T];L^nR)) < CR-"^, 

(7) ||V(^x - ^.«)||i.o([o,T];L2(Q)) <Ci?-/5. 
Above, a and [5 are given by: 

{ 1/2-0/2 i/m/o' ^^-^^ 
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and 

( 1/2 + 30/2 ifm = 0, 
^~ \ 1/2 + 0/2 ifm^O, ^'^■'> 

in two dimensions, while a = 1/2 and (3 = 3/2 in three dimensions. 

For case I of Definition the constants above depend only on 0,; for 
cases II and III some of the constants also depend on T . 

4. Energy argument 

Proof of Theorem The proof proceeds much as in Section 2 of [5] : 
Using our approximate solution to {E) we make an energy argument to 
bound the difference 

W = u"'^ - 

in the norm. Then using inequality (6) of Proposition 13.11 we apply 
the triangle inequality to complete the proof. We give the argument in 2 
dimensions only; it is valid with minor adaptations in 3 dimensions. The 
only delicate point in adapting to 3 dimensions is that we deal with weak 
Leray solutions for which we cannot perform energy estimates. However the 
energy inequality is equivalent to the necessary estimates; see |5j for a more 
detailed discussion of this issue. 

Subtracting the identity in (j3.4p from (NS) we obtain 

dtW - i^AW = -u^^^-Vu"^^ - Vp"'^ + uAu^ + ^p^u ■ Vu 

+ ip^Vp - dttpV-^ip^'. 

Multiplying both sides by W and integrating over Qji gives 

\-^\\W\\l2+u\\VW\\l2 =Ii + h + h + h + h, 
where when s > 1, 

^3 = Ln^^'^u- Vu) -W, h = f^^ i^Vp ■ W, 

In Ii we integrated by parts to remove Au^, but when s > 2 it is more 
advantageous to retain it, using 

h = u f W ■ Au^. 

When s > 1 we apply the Cauchy-Schwarz and Young's inequalities to 
the first form of Ii to get 

< ^ (||VW^||^2 + llVn^ll^^) < ^ ||VT^||^2 + Cv, 
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and when s > 2 we apply the Cauchy-Schwarz inequahty to the second form 
of Ii to get 

Summing I2 and I3 and using 

[ (u^'^ • VT^) ■W = 

we have 
I/2 + /3I 

{u^ ■ V(n - u^)) -W -{W ■ Vu^) ■ W + [{if^u - u^) ■Vu]-W 



/< 



< Wu^'Wl^ p{u - n«)||^, lliyil^, + IIVu^IIl- \\w\\l, 

+ y^U-U^L2\\Vu\\L^ \\W\\j^2 

< CR-P \\W\\^2 + C WWWl^ + CR-'' \\W\\l2 

<c{R-<' + \\w\y)\\w\y, 

where we used inequahties (6) and (7) from Proposition 13.11 and also that 
a<p. 

Summing I4 and I5 and integrating the first term by parts gives, 



\h + h 



^R 



< (iIpV^'^IIl^ + imv^v^'h^) \\w\\^, < CR-' \\w\\^, 

<CR-^^ + C\\W\\l,. 
When s > 1, we conclude that 
ld_ 
2di' 



-Mw\\l, + u\\vw\\l, 



< ^ \NW\\l2 +Cu + CR-" \\W\\l2 + CR-^^ + C \\W\\l2 



Integrating in time gives 



< \\W{0)\\l2 + CTv + CTiR-^"" + R-"^^) + C [ 

Jo 



2 

L2 
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It follows from Gronwall's lemma that 

II^IIlo=([o,t];L^(q«)) < {mf + C{u + R-^^ + R-'')) ^'^ e^^ 
< {f{R) + C{v^'^ + + R-')) e^^. 
Then from the triangle inequality and inequality (6) of Proposition I3.H 

II V R II 

\\U ' - u\\L°°ilO,T];L^(nR)) 



< \\u 



R 



^llL°°([0,T];L2(nfl)) + ll^llL°°([0,T];L2(nfl)) 



< ( F{R) + C(z.i/2 + R- + R-')) e^^. 



When s > 2, we have instead that 

< Cu \\W\\l2 + CiR-"" + R-^) \\W\\j^2 + C \\W\\l2 . 

Dividing both sides by ||W^||l2 (it is easy to see that division by zero will 
not invalidate the following inequality after integrating in time) gives 

^\\W\\l2 < Ci^ + CR-"" + CR-^ + C \\W\\r2 . 
Integrating in time and applying Gronwall's lemma, we have 
||1^||l2 < (f{R) + C{u + CR-"" + R-^)] e^^. 



The bound on \\u^'^ — if||Loo([o^T];L2(f7fl.)) follows from the triangle inequality 
as for s > 1. 

The value of a in the statement of Theorem ll.2l is chosen so that a = a{6) 
gives the optimal rate of convergence in each case; this corresponds to = 1 
for Case I; 9 = 1/3 for Case II so that 6 = a; and 9 = 1 was fixed for Case 
III. □ 

5. Truncation operator in 2D 

Let u be in Em H for some m in M with vorticity u!{u) having compact 
support in a ball of radius Rq. Let ^jJ be the stream function, as defined by 
the expression in (j2.ip . 

Let ip^ and be defined as in Section [3] and recall the definition of 
and Tr given in (j3.2p . = Tru = S/-^{ip^il)). To explore the properties of 
Tr we must first establish some bounds on the norms of u, Vu, and if) in 
S/j. To this end we write 

u = V + a, 

where v is in Vc{^) and a is a stationary solution with radially symmetric, 
smooth, compactly supported vorticity of integral m; we assume that w(cr) 
is of distinguished sign. We can assume, without loss of generality, that 
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the support of w(ct) is also contained in the ball of radius Rq, from which 
it follows that the support of uj{v) is contained in this same ball. Now, 
V and (T are also divergence-free vector fields and hence we can define 
their associated stream functions and ipa using the expression in (j2.1l) . 
But then v = S/'^ipy and a = V'^ifj^- It follows in particular that v can 
be written in terms of uj{v) through the Biot-Savart law v = K * an 
integral operator with kernel 

From this explicit expression and using the fact that the integral of uj{v) 
vanishes, together with the easily obtained estimate 

< 2||a;(u)||Li(M2), 

it follows that there exists C = C[Rq) > such that 

\v{x)\ < CCo/ \xf , |V^;(x)| < CCo/ \xf (5.2) 

for all |x| > 2i?o, with Co = 2\\uj{u)\\ii^^2y Similarly, it follows from the 
explicit expression for ^py, (j2.ip . that 

\tpv{x)\ < CCo/ \x\ for all \x\ > 2Ro. (5.3) 

Put coordinates on S/j as in the definition of (pn in Section [3l Letting a 
be the length of Fi it follows that the length of F/j is aR. Then 

UvWh^j^^) = \J{s,r)\ \ij^{s,r)\^ drds, 

where J is the Jacobian of the transformation from rectangular coordinates 
to {s, r)-coordinates. Because of the way we constructed S/j and because 
1^ ^, \ J\ < C and S/j lies outside a ball of radius C{Q)R. Thus by (j5.3p . 
\'ipvis,r)\ < CCo/R in the integral above as long as C{Q)R > 2Ro; that is, 
as long as 

R > fi{n)Ro, (5.4) 
where /x(r2) = 2/(7(0) depends only on the geometry of 17. Then 

UvWl^j:^) < (cClR-^aRSiR'Y^ < CCoR'^^-'/^. 
Since each of V"*" and V introduces an extra factor of 1/ \x\, it follows that 

< C7Coi^^/^-3/^ l|Vr.||x.(E«) < CCoR'/'-'/'. 
As previously pointed out, see (jl.ip . a is given by 

/ X2 Z"''^' xi fl^'l \ 

a{xi,X2) = -—2 / ruj{a){r) dr, / ruj{a){r)dr 
\ \x\ Jo \x\ Jo I 
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SO that \cr{x)\ = \m\ (27r) ^ / \x\ for jx| > Rq. Thus, a decays hke ipy so we 
can see that 

Mm^^) < C \m\ ||Va||i2(£^) < C \m\ R'/^-^/^. 

The expression for ip^- can be calculated directly using (j2.ip together with 
the radial symmetry of oj{a). Of course, we can add an arbitrary constant to 
ipa and still satisfy the equations a = V-^xpr^ and Aipf^ = u!{a). For > Rq 
we obtain: 

V^ix) = — log|x| + C. 
zvr 

Since when m 7^ we assume that is a disk centered at the origin, 
we can choose the constant Cr so that ipf^ = on Tr. The value of dt^ 
is unaffected by the choice of Cr, however, and dttp in inequality (4a) of 
Proposition 13.11 is the only direct use of tp that we make, so the choice of 
Cr, though it depends on R, will not affect any of our estimates. 

Applying Poincare's inequality (or integrating ■0o- directly) gives 

The factor of R^ here comes from the thickness of Tir. 

Adding the corresponding bounds for u = v + a and ip = ijjy + ipa, 

UhH^n) < CCoR'/'-'/' + C \m\ R^'l^-^l\ (5.5) 
and for the velocity 

II^I1l2(e«) < CC^R'!''-''!'' + C \m\ R'I^-^I^ (5.6) 

and 

l|Vn||,.2(s^) < CCoi?^/2-^/2 + C \m\ r"^-^'\ (5.7) 

These inequalities each hold as long as (15. 4p holds. 

Let X be the subspace of all vector fields in Em ni?i(M^) whose vorticity 
has compact support. We can now describe the relevant properties of the 
two dimensional truncation operator, adapting Lemma 4.2 of [8j. 

Proposition 5.1. Let Q be a disk centered at the origin and let the trun- 
cation operator Tr he defined as in 113. 2\) . Then Tr: X ^ V{i^R) with the 
property that for all u in X 

\\u - Tru\\l2(^^) < C||a;(n)|iii(K2)i2~" (5.8) 

and 

\\V{u - TRu)h2inn) < CMu^mu^^R-" (5.9) 

for all R satisfying (5^, where C = C{^}), a is defined in (3. 6\) . and jS is 
defined in \3. Tj ). Also, 

\\^Tru\\l^(^r) < C\\uj{u)\\ 12(^^2), (5.10) 

where C = C{Rq). 
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Proof. Using the inequahties in ()3.ip . (|5.5p . (15. 6p . and (15. 7p . we have, 

< 11(1 - <^^)V^^||i2(^^) + m^if^'h^inn) 

< IkllL^cs^) + I|V(^^|Uoo(s^)||V'||l2(e^) 

< CCoi?''/2-3/2 + C \m\ i?^/2-i/2 + C'ii-^Coi?^/2-i/2 

+ CR-' \m\ i23e/2-i/2 

+ C\m\R'/^-'/^ 
<CCoR-"^-^'^+C\m\R'l^-^'\ 

In the final inequahty we needed to only include the two terms that dominate 
(depending on whether m = 0) for all in [0, 1]. Clearly, \m\ < \\ijj(u)\\ii(^2-j, 
so this gives (|5.8p . 
Similarly, 

llVn - VTRu\\L2ina) = " VV^(99«V')||i2(f,^) 
= \\Vu - V(v9^V^V) - V(V'V^(^«)||i2(f,^) 

= iivti - ip^vv^i! - Vip^ (g> v^i) - vv- «) v^ip^ - V'VV-^v3^||i2(f^^) 

= 11(1 - if^)Vu - Vip^ V^V - VV' O V^(^^ - V'VV^(/?^||i2(n^) 

< I|Vu||l2(2^) + 2 ||VV9^||^^(^^^ ll^i||L2{E«) + l|VVv9^||i^(s^) ||V'IIl2(2^) 

< CCoR'^'-'/' + C \m\ r"^~^I^ + Ci?-^ (Coi?^/2-3/2 + \m\ R'/^'^'^) 

+ CR-^' [CoR'^^-'/^ + \m\ i?3^/2-i/2 j 

= CCoi?^/2-^/2 + C |m| i?^/2-3/2 + C'Coi?-''/2-3/2 + ^ ^-e/2-1/2 

+ CCoR'''/'-'/' + C\m\R-'/^-'/' 

< CCoR-^'/^-'/^ + C \m\ R-'/^-'/\ 

Since, as we noted, \m\ < \\u){u)\\ii(^2-^, the estimate above gives (|5.9p . 
The bound on || VTRn||^2(f2^) is obtained in the same way except that 
l|Vn||i2(s^) is replaced by || V'u||i2(Q^\s_R) < l|Vn||i2(R2) = \\u{u)\\L2^^2y 
But ||w||li(r2) < C||a;||i2(-iR2) because the support of w is contained in -B/Jq, 
giving ([5TU]) . □ 

Remark 5.2. Had we not added the restriction that 17 is a disk, then since 
ipa- is defined only up to an additive constant, we could choose this constant 
so that f^^ Vo- = 0. By normalizing tp^j this way, though, the factor of R^ 
in Poincare's inequality becomes the diameter CR oi VLr and we have 

IIV'.IIl^ce^) < Ci?|k||i2(2,) < Cmi?^/2+i/2. 
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This gives 

UWLH^a) < CCoR'"-'" + Cmi?^/2+i/2 (5.11) 
and we would have had the following bound in the proof of Proposition 15.11 

< CC7oi?^/2-3/2 ^ ^^^e/2-1/2 ^ cCoR-"^-^'^ + Cmi?-^/2+i/2 

<CCoR"^-^'^ + CmR-'l^+^l\ 
No value of 6 in [0, 1] will allow this to vanish as i? ^ oo when m ^ 0. 

6. Estimates in 2D 

Proof of Proposition \3.1\ in 2D. Define n^, ip^, and S/j as in Section[3l 
Then 

= V^{^^^) = ^^u + ^V^^^ (6.1) 

where ■0 = V'l' + V'o- as in Section [5l 

The and norms (indeed, all norms) of the vorticity for solutions 
to {E) are conserved over time, while the bounds on the (for m = 0) 
and norms of the velocity are bounded over any finite time interval 
(the former is conserved, and a bound on the latter is uniform). Thus, 
the estimates in Proposition 15.11 when applied to our solution u to {E) are 
uniform in time. Hence, in the estimates that follow we will generally not 
explicitly refer to the bounds in time. 

Inequality 1: follows directly from Proposition 15.11 

Inequality 2: We have 



R\ 



From (|5.3p . llV'f IIloo(Sjj) ^ CCq/R and applying Poincare's inequality in the 
L°° norm, 

UJi^i^n) < C'«'lkllL-(E«) <C\m\R'R-' < CCoR'-\ 
We conclude that inequality (2) holds for sufficiently large R. 



Inequality 3: We have, 

Vn||i«>(n^) + 2 



Vu-^IIloc < ||Vn||i^(o^) + 2||V(p^||^^(^^) ||n||ioc(s^) 



+ l|VV if'^h^ij^j,) IIV'IIl-(Sb) 
<C + CR-'^ + CR-'^'^CqR'^-^. 
We conclude that inequality (3) holds for sufficiently large R. 
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Inequality 4a: We begin with the observation that 

A-ipt = -w(div (n (g) n)) (6.2) 

and 

— Ap = div div (u ® u) . (6-3) 
Now, u is bounded, uniformly over a finite time interval, in L'^{M?). To see 
this write u = v + a; clearly a G L^(M2) and v G L°° n L2(M^). It follows 
that the right-hand side of both (|6.2p and (j6.3|) are second derivatives of an 
vector field. We can express Tpt and p as linear combinations of Riesz 
transforms of terms which are uniformly bounded in L^(]R^) and hence, by 
the C alder on- Zygmund inequality, it follows that ipt and p are both bounded, 
uniformly over a finite time interval, in L^(IR^). Thus, 

< \\p\\L^j:n)\\^'P^h^{^R) + ll^t^llL2(Efl)llVi/3-^||L-(s«) 

< CR-^ + CR-^ = CR-^. 

Inequality 5: For / a scalar and v a vector field we have that 
A{fv) = Afv + 2V/ • {Vvf + fAv, 

so 

Au^ = A{ip^u) + A(V'V^(/?^) 

= Aif^u + 2V(^-^ • {Vuf + (^-^Ati + AV'V-^(^^ 
+ 2VV' • (VV^(^-^)^ + ^AV^99^. 
Then because Ai/' = a; and u = S/'^tp, 

< C ||An||^2(iR2) + C, 
for sufficiently large R by the bounds in (15. 5p . (j5.6p . and (15.71) . 

Inequality 6: We have, - n||i2(s^^) < CCqR'"' by dSSl) and 
ll-u^ - V9^n||i2(f^^) < - n||i2(f^^) + 11(1 - (/?ij)'u||i2(f^^) 

< ||n^-n||i2(n^) + ||u||i2(s^) 

< CCoR''" + CCoi?^/^"^/^ + C [m| i?^/2-i/2 

< CCoR^'' 
where we used (15.61) and dH 
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Inequality 7: the two-dimensional case follows directly from Proposition l5.ll 

□ 

7. Decay of velocity in 3D 

The basic existence result for a solution to {E) is given in Theorem 17.11 

Theorem 7.1. Assume that uq is in C"^(IR^) n Vc'(IR'^) for s > 1. Then there 
exists T* in (0, oo] such that for all T in (0, T*) the solution u to (E) in the 
whole space lies in L°°([0, T]; C^(M3)) yjUh Vp in L°°([0, T]; L'^(IR3)) for all 
a in (l,oo]. Also, the vorticity oj = uj{u) lies in L°°([0, T]; L"(]R'^)) for all a 
in [l,oo] with a hound on its norm that is independent of a. Furthermore, 
IV remains compactly supported for all time, with the support contained in 
a ball of radius R{T) < + ||'"I1loo(-[q ^-JxR^) ^ /o*" '^^^ ^ [0;^]; where the 
support of loq is contained in Bfi^^{0). 

Proof. The initial vorticity ui = uj{uq) is compactly supported and continu- 
ous by assumption so Vfo is in L"(M'^) for all a in [1, oo]. Thus by Theorem 
4.2.3 p. 79 of d], there exists T* in (0,oo] such that for all T in (0,r*) the 
solution u to (E) lies in L°°([0, T]; C"(M3)) with Vp in L°°([0, T]; L"(M3)) 
for all a in (1, oo). 

Taking the vorticity of {E) gives 

-j^ = OtuJ + u ■ vuj = uj ■ Vn, 

where Dco/Dt is the rate of change of the vorticity along the flow lines. 
Assume that the support of a;'' is contained in a ball -Bij(o) of radius R{0). 
Then until some flow line starting from within -Bij(o) reaches a point outside 
of -B^(o) there can be no change in vorticity at that point. This gives the 
compact support of the vorticity for all time along with the bound on its 
support. □ 

Remark 7.2. Theorem 17.11 is stated for three dimensions but holds in two 
dimensions as well, with T* = oo, with minor adaptations to allow for 

uo G C'{R'^)nEm, m / 0. 

Corollary 7.3. Let u be a solution to (E) as in Theorem \7.1\ Then u is in 
L°°([0,r];LP(R3)) for all p in (3/2, oo] and Vu is in L'^{[0,T]; LP{R^)) for 
all p in (1, oo]. Also, 

\u{x)\ <C/\x\'^ , |Vu(x)| < C/|x|^ 
for all t in [0,T] and \x\ > 2R{T), and 

Ht)\\LH^n) ^ CR-'/^ l|Vn(t)||^.(2^) < CR-'/', 
for all t in [0,T] and all R sufficiently large that S/j is contained in B^j^^rpy 

Proof. Let \I' be the associated stream function in three dimensions and 
recall the expression for ^ given in (j2.2p . In particular, u = curl ^ and 
— A^f = (jj, where lo = uj{u) is the associated vorticity. This gives rise to 
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the three-dimensional Biot-Savart law, allowing us to write the velocity in 
terms of vorticity as: 

— 1 f X — y 

u = u{x) = — r^xuj{y)dy. (7.1) 

Air J \x — 

The decay of u and Vii then follow from the compact support of to and the 
Biot-Savart law. Because u and Vu are in L°°([0,T] x W^), the membership 
of u and Vu in the stated spaces and the decay of their norms then follow 
from their decay at infinity. □ 

8. Truncation operator in 3D 

Let ip^ and S/j be as in Section [3] and recall the definition of the truncation 
operator in three dimensions: if u in H^{R^) then Tru is given by (jS.Sp : that 
is, Tru = V X [if^"^). Set E = E{x) = (47r|x|)~"^, the fundamental solution 
of -A in R3. 

Lemma 8.1. // / in L^(]R^) n L°°(R'^) is compactly supported in the hall 
Bl of radius L centered at the origin then for all R> R* , 

\\dkE*f\y^j,^^<CLy'\\f\y^^,^R~'/\ 

ll^*/llL=^(E,)<C^^'/'ll/llL^(M3)i^~^ 

where C is an absolute constant and R* depends only upon Qi and L. 

Proof. On the compact support of / gives \E * f\ < C ||/|lii(K3) / \x\ < 
CL^/'^ ll/llL2(ig3) / \x\. Therefore, as long as R is large enough so that is 
wholly contained in i?^, 

\\E* f\\L2^Ba)<CL''^\\f\\L^(M:^)R''\ 
and similarly for the L°° norm and for the inequalities for dkE * f. □ 

Proposition 8.2. Let u be a solution to (E) as in Theorem \7.1\ Then there 
exists R* > such that for all R> R* and all t in [0,T], Tru is in V{Qr), 
and 

\\u-Tnu\\L2(^n^)<CiR-^/^, (8.1) 

||V(u-THn)||^2(f,^) <C2i?-3/2. (8.2) 
The constants Ci and C2 depend on u and T, and R* depends on RiT) and 
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Proof. On Tji, T^u vanishes identically, while clearly div Tnu = on Qjj. 
Now, Tru = ipRU + VipR X so 

i,k 

We used Corollary 17.31 and Lemma 18.11 in the final inequality. 
Also, 

dkTRU = (fRdkU + VdkipR X ^ + {dk'PR)u + V(pR x dk"^. (8.3) 

Hence, 

\\VTru - Vn||i2(n^) 

< l|Vn||z.2(s,) + II VV(^«||^^(^^) miLH^^) 

+ I|vv'''ILoo(s,)Eii^^^*^(")IIl^{2«) 

A: 

A: 

where we used Corollary 17.31 and Lemma l8.ll in the final inequality. Thus, 
Tru e V{i^R), as desired. □ 

9. Estimates in 3D 

Proof of Proposition \3.1\ in 3D. Define ip^ and Tjr as in Section[3l Then 

= Tru = ipRU + VipRX (9.1) 

Inequality 1: follows from ()8.2p . 
Inequality 2: From (j9.1|) . 

^ II^IIl«=(r2) + Ci^'-^IIE; * w||l°°(Sb)' 
which, using Lemma |8.H is enough to establish inequality (2). 
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Inequality 3: Recall the expression for Vu-^ in ([83]). Then, using Lemma FS-H 
we have 

llVn-^llioo < ||V'u|[ioo(Qii) + |iVV99-^||ioo(2fl) |j*|lLcx,(sfl) 
+ II'^V'^IIl°°(s«) W^kE * u;||^oo(sfl) 

k 

<C + CR-^ + CR-^ + CR-^ < C. 

Inequality 4b: The proof of this inequality is similar to that of inequality 
4a. We begin by observing that satisfies an equation analogous to (|6.2p . 
namely 

A^'f = curl div u®u. 
Now, u € L^(]R^) n L°°(IR^), uniformly in time up to T* . Therefore, us- 
ing again the Calderon-Zygmund inequality, it follows that ^'t G L^(]R^), 
uniformly in time up to T* . The estimate for the pressure p in three dimen- 
sions is exactly the same as for two dimensions, using (I6.3p . Thus, using the 
scaling of (pR, we get 

L2(n_R) 

< l|p||L2(SB)l|Vv'^||L-{Efl) + l|5t^||L2(s^)||V(^-^||i^(s^) 

< CR-^ + CR-^ = CR-\ 

Inequality 5: We must compute An^; for the sake of simplicity we keep 
track mostly of the order of derivatives. We have then that 

An^ = {A(p^)u + 2Vip^ ■ Du + f^Au + Vip^ x to 

+ 2D'^ip^ ■ D-i! + VA(p^ X ^. 

Using Corollary 17.31 and Lemma 18. H we find 

+ C ||A'u||^2(f^^) + CR-^ II^IIl2(e^) 

+ CR-'\\DnLHJ,,) + CR--'\mLH^,) 
<C + CR'^/^ < C. 

Inequality 6: We have, \\u^ - u\\l2(^^^^ < CR'^/"^ by dSTT]) and 

\\u^ - V?^n||i2(f^^) < \\u^ - n||i2(f^^) + 11(1 - ipR)u\\L2^^^) 
< \\u^ - u\\l2^^^) + |ln||i2(2^) < CR~^'^, 
where we used Corollary I7.3[ 
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Inequality 7: the three-dimensional case follows directly from Proposi- 
tion [821 □ 

10. Truncation of the initial velocity 

The following observations hold in any dimension. 

Define Wr: F(M'^) H{nR) as follows: for any u in y(R'^) let Wru be 
that unique vector in H{i}R) having the same vorticity as u (ft is simply 
connected so the vector is unique) , see the paragraph following the statement 
of Theorem II .21 Let -Ph-(q^) : y(M^) H(i^R) be the projection operator 
into H{nR). 

Lemma 10.1. The operators Wr and PH{nji) o-Te the same. 

Proof. Let u be in ^(M'^). Then by the Leray-Helmholtz decomposition 
there exists v in H{Qr) n -ff^(r2/j) and p in such that u = v + Vp 

on ^}r, with Ap = and with v and Vp unique. Then v is orthogonal to 
\7p in L'^(Qr) (that is, in the H(QR)-noicm) so v = Ph'(q^)U. But the curl 
of a gradient is zero, so uj{v) = uj{u) which shows that v = WrU as well. □ 

Corollary 110.21 follows immediately from Lemma 110. H given the well- 
known properties of the projection operator. 

Corollary 10.2. Wr has the same convergence properties in the LP'{Q.r)- 
norm of the velocity as does Tr. In particular, we obtain (II. 2p . 

Remark 10.3. The operator Wr is a more natural "truncation" operator 
than Tr since it involves literally truncating the vorticity. It is perfectly 
adequate for truncating the initial velocity for the Navier-Stokes equations, 
but because we do not know how to control the remainder following trunca- 
tion in the /7^-norm, we cannot use it to construct the approximate solution 
to the Euler equations. 

11. Comments and conclusions 

As a first comment, we note that the estimates in Sections [7] and [8] for 
three dimensions are considerably simplified by the assumption that the 
initial vorticity is compactly supported, whereas for two dimensions this 
assumption is merely a minor convenience that allows us in Theorem 11.21 
to give an explicit rate of convergence in R. If one drops the assumption 
of compact support, it does not seem possible to obtain a uniform-in-time 
bound on the decay of the gradient of the velocity. One can obtain a bound 
on the decay of vorticity, however, if one assumes that the norm of the 
initial vorticity on a ball of radius R decays for large R at least as fast as 
CR~^/^. One then modifies the energy argument in Section HI integrating 
by parts differently to use the decay of the vorticity in place of the decay of 
the velocity. The value of a in Theorem 11.21 however, must be strictly less 
than 1/2. 
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The case of initial vorticity with nonzero integral in two dimensions corre- 
sponds to the situation where the limiting full plane flow has infinite energy. 
Since our argument is based on energy estimates, it is natural that this situ- 
ation would be complicated. In this situation we have only studied the case 
of the expanding disk, but, as explained in Section 5, this is not a matter 
of convenience. Our argument makes essential use of the fact that we are 
working in a disk. However, the restriction to the expanding disk is by no 
means natural, and removing this assumption becomes an interesting open 
problem. 

Finally, we note that our result was proved for smooth flows but that, in 
two dimensions, it would be very reasonable to consider initial vorticities in 
L^, p > 2, compactly supported, given that, in our proof, the high regularity 
was needed in estimates near the boundary, where vorticity acts as far field. 
If p > 2 then the support of vorticity can be controlled since the Euler 
velocity is a priori bounded. 
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